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The penetration depth is calculated over the entire doping range of the cuprate phase diagram 
with emphasis on the underdoped regime. Pseudogap formation on approaching the Mott transition, 
for doping below a quantum critical point, is described within a model based on the resonating 
valence bond spin liquid which provides an ansatz for the coherent piece of the Green's function. 
Fermi surface reconstruction, which is an essential element of the model, has a strong effect on the 
superfiuid density at T = producing a sharp drop in magnitude, but does not change the slope of 
the linear low temperature variation. Comparison with recent data on Bi-based cuprates provides 
validation of the theory and shows that the effects of correlations, captured by Gutzwiller factors, 
are essential for a qualitative understanding of the data. We find that the Ferrell-Glover-Tinkham 
sum rule still holds and we compare our results with those for the Fermi arc and the nodal liquid 
models. 
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I. INTRODUCTION 

The initial discovery of superconductivity in the 
cuprates precipitated a rush to find higher values of the 
critical temperature Tc. In this onslaught, not only were 
new superconducting members of the cuprate family dis- 
covered but it was quickly realized that, through oxygen 
doping or doping with other elements, a particular com- 
pound could display a range of Tc values. With hole- 
doping, in particular, these Tc values could be quite high 
and showed a dome for Tc as a function of doping in 
the phase diagram. While studying optimal doping for 
the maximum Tc became the primary focus in initial re- 
search, it was later appreciated that the unusual phase 
diagram of these materials was of great interest in itself. 
Understanding why a maximum Tc exists and what con- 
trols the reduction of Tc away from maximum is hoped to 
elucidate the physics of the interactions involved. Ideally 
this would give direction for what might be attempted in 
materials development in order to enhance Tc- 

Further research into the general phase diagram oc- 
curred on many fronts, with evidence of an antiferromag- 
netic (AFM) insulating state in the parent compound and 
at low doping, and strange metallic behavior in the nor- 
mal state above the superconducting dome. The most 
intriguing discovery is possibly the presence of what is 
now termed a "pseudogap" feature, occurring in the nor- 
mal state on the underdoped side of the superconducting 
dome, which exists for higher temperatures as the AFM 
state is approached.^ This pseudogap is an energy gap- 
like feature seen in normal state properties which has 
led to a number of imaginative theoretical proposals for 
its existence. Superconductivity in the cuprates is now 
thought to have d-wave superconducting pairing, most 
likely due to spin-fluctuations."3'»2, It is thought, there- 



fore, that the approach towards the AFM state should 
enhance the spin-fluctuation pairing interaction V and 
hence increase Tc. However, if this pseudogap represents 
a competing phase, possibly associated with the AFM- 
Mott insulator, it could be responsible for a reduction in 
the electronic density of states at the Fermi level N{0) 
which could then reduce the superconducting Tc which 
depends on N{0)V in simple BCS theory. Consequently, 
it is natural that many proposals for the pseudogap state 
have been based on a competing phase, such as, the d- 
density wave theory.** Others have been suggested which 
include the ideas of preformed pairs^iii^ arising above 
T 

The pseudogap phase and how it may affect supercon- 
ductivity has become a major focus of both theoretical 
and experimental work. Even knowledge of where the 
pseudogap line in the phase diagram might end [possibly 
at T = and possibly as a quantum critical point (QCP)] 
is still a matter of debate. Some experiments suggest it 
ends on the edge of the superconducting dome on the 
overdoped side and others place it ending at T = inside 
the superconducting zone anywhere from optimal doping 
to overdoped.^ It should be noted that other works on 
QCPs in heavy fermion superconductors usually suggest 
such a point existing under the superconducting domeJ^ 

In spite of mounting experimental work on the under- 
doped cuprate superconductors, it has been hard to de- 
velop a microscopic theory that can build in the physics 
associated with the approach from the metallic state to 
the Mott insulator and give a workable formalism, which 
includes doping, for providing theoretical insights and 
predictions for experiments both in the normal and su- 
perconducting state (including the idea of a pseudogap) . 
Further to this issue has been the more recent experimen- 
tal development where a possible reconstruction of the 
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Fermi surface in the underdoped cuprates has been found 
to occur based on seeing Fermi pockets from de Haas- 
van Alphen experiments^i and or Fermi arcs in angle- 
resolved photoemission (ARPES) experiments.^^ There 
has also been one report of an observation of pockets 
from ARPES.— Moreover, there has been indirect ev- 
idence of arcs from other experiments, such as optical 
spectroscopyi^, specific heat^i^ and scanning tunneling 
spectroscopjii^ measurements. Indeed, attempting to re- 
solve the arc versus pocket debate has been the subject 
of numerous papers. 

One theory which is promising in this regard has been 
due to Yang, Rice and Zhang— (YRZ) which is based 
on previous studies of a resonating valence bond (RVB) 
spin liquid state originally proposed by Anderson^^. The 
merit of the YRZ approach is that they have used pre- 
vious numerical and theoretical RVB-type studies to de- 
velop an ansatz for the many-body Green's function that 
would represent the RVB spin liquid. This ansatz builds 
in the approach to the AFM-Mott insulator and shows a 
reconstructed Fermi surface which is a large Fermi sur- 
face when there is no pseudogap, above optimal doping, 
but is reconstructed to small Fermi or Luttinger pockets 
(which look more like arcs when the quasiparticle weight 
is included) for the underdoped case. The pseudogap 
opens up around the AFM Brillouin zone in this theory 
and it gaps out or reconstructs part of the Fermi surface. 
So far this theory has been used to evaluate a number 
of experimental propertie a^^'^°i^^'^^ with good qualita- 
tive agreement and, in the same spirit as BCS theory, 
using this ansatz allows us to determine what are the es- 
sential elements that should go into a more sophisticated 
microscopic theory should it be developed in the future. 

In this work, we wish to examine the long-standing 
puzzle associated with the penetration depth measure- 
ments. The penetration depth was one of the first ex- 
periments to clearly indicate that a d-wave order pa- 
rameter symmetry for the superconductivity was present 
in the cuprates This was an extremely important re- 
sult in influencing the direction of research in this field 
as it allowed for the elimination of a number of possi- 
ble mechanisms for Cooper pairing. It also clarified the 
need for very high quality samples to remove the obscur- 
ing features due to impurities. The essential observa- 
tion from the experiment was that the superfluid density 
Ps{T), which is related to the penetration depth A(T) by 
PsiT) oc 1/A^(T), showed a low temperature linear T be- 
havior as expected for a clean BCS d-wave superconduc- 
tor, i.e., Ps{T) = Ps(0) — bT. However, in the underdoped 
regime, it is known^iiS^i^^ that while the zero tempera- 
ture value of the superfluid density depends strongly on 
the doping x, the coefficient b of the first linear-in-T cor- 
rection is much less sensitive to x. This result cannot 
be understood within a simple BCS d-wave model. Our 
goal in this paper is to study the penetration depth in 
the YRZ model and to see if the experimental data, with 
its doping dependence, can be explained by this theory. 
Furthermore, we wish to see if there is evidence for a re- 



constructed Fermi surface in the penetration depth data. 
With this study we can develop a better understanding of 
the physics which is giving rise to this non-BCS behavior 
of the doping dependence of the penetration depth. 

Our paper is structured as follows. In Section II, we 
introduce the basic features of the YRZ theory that en- 
ter into our calculations. This is followed by a discus- 
sion of the penetration depth formula used in this work 
and its various limits, given in Section HI. In Section IV, 
we summarize some theoretical formulas associated with 
the frequency-dependent optical conductivity, a quantity 
which is also related to the penetration depth. This will 
aid in our discussion of the Drude weight and the sum 
rule. We then present our results in Section V and pro- 
vide our conclusions in Section VI. 



II. THEORETICAL MODEL OF YRZ 

The YRZ model provides an ansatz for the coherent 
part of the many-body Green's function for the case of 
a doped RVB spin liquid. It includes a dependence on 
doping X and is given asj^i^i 

where 

TTilt _|_ TP 

- 2 

~ i^k+ek) 

t;k - 2 ' 
»^ = i(l±|). (2) 

In these expressions ^2 — —2t{x)(coskxa + coskyo) 
is the first nearest neighbor tight-binding dispersion, 
which introduces the antiferromagnetic Brillouin zone 
boundary (AFBZ) into the Green's function along which 
the pseudogap Apg develops. The band structure 
^fe = —2t{x) {cos kxtt + cos ky a) ~ -it' (x) cos kxd cos ky a — 
2t" {x){cos2kxa + cos2kya) — fip is taken from tight- 
binding for a system which includes hopping terms up 
to third nearest neighbor, with ^p, a chemical poten- 
tial determined by the Luttinger sum rule.^^ Doping de- 
pendence enters in three ways into the Green's function. 
First, the band structure is doping dependent through 
the hopping coefficients: t{x) ~ gt{x)tQ + igs{x)Jx/S', 
t'{x) = gt{x)t'Q, and t" [x) = gt{x)t'^, where gt{x) = 
2x/ {l + x) and gs{x) — 4/(1 + a;)^ are the Gutzwiller 
factors and J/to = 1/3 and x = 0.338. YRZ use tg/^o = 
—0.3, t'^/to ~ 0.2, a choice of parameters to match 
this energy dispersion to that calculated for Ca2Cu02Cl2 
(Ref.^^). The x-dependence of these coefficients reflects 
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the fact that strong correlations will narrow the bands as 
the Mott insulator is approached. Second, the coherent 
part of the Green's function changes with doping with a 
Gutzwiller factor given by gt = gt{x). Such factors re- 
flect the projection out of doubly occupied states and the 
approach to the atomic limit. Third, the magnitude of 
the pseudogap, Apg, and superconducting gap. Age, are 
also doping dependent, as inferred from experiment and 
the behavior of Tc. These gaps are taken to be d-wave, 

A° (x) 

such that we can write them as Apg = — e| — [cos k^a — 

coskyu) and Age — ■^^^^(cosfcj^a— cosfcyo), respectively, 
with a the lattice constant, Apg(x)/io = 0.6(1 — x/0.2) 
and A°^{x)/tQ = 0.14(1 -82.6(a;- 0.2)2). general, the 
gap could contain many higher harmonics^ as is also 
the case in conventional superconductorsSSiSSii^S but such 
complications are not essential for a first understanding. 
From Eq. ([l]), one can extract the YRZ spectral func- 
tion A{k, Lu) and see that there are four energy branches, 
given by the energies, ±i?5, where Eg = y^iE^)^ + A^^. 
Following the usual development of the equations of su- 
perconductivity, both the regular and anomalous spectral 
fmictions arc found to be: 

A{k,cu) = J2gtW^[{u-f6ic.-E^) + {v"fSic. + E^^)], 

a=± 

(3) 

B{k,Lu) - 9t^>"^"i^('' - ^s) - S{co + E^)], (4) 



fer to a Fermi liquid which is just the case of taking the 
pseudogap to be zero in our formalism. We also use the 
term YRZ modified (YRZ mod.). In the YRZ model, the 
presence of the pseudogap reconstructs the large Fermi 
surface (found on the overdoped side of the phase dia- 
gram) into a small Fermi pocket. If we use the YRZ 
model but take the superconducting d-wave gap to be 
finite only in the region of the Fermi pocket and to be 
zero beyond the pocket towards the antinodal direction, 
then we will call this YRZ modified. In conventional 
superconductors^SiSSi^ such as Pb and Al, there are di- 
rections where the Fermi surface is gapped out by the 
crystal potential and one finds that the superconduct- 
ing gap is also zero there. In the present case it is the 
pseudogap which prevents the superconducting gap from 
having its full amplitude in certain regions. 

Two other models in the literature are the nodal liquid 
and the Fermi arc model. In these two models, the pseu- 
dogap is taken to be on the large Fermi surface given by 
^fc. In the nodal liquid, both the pseudogap and super- 
conducting gaps are active over the entire Fermi surface. 
In the arc model, the pseudogap is only finite in a region 
of the Fermi surface near that antinodal direction on an 
arc of the large Fermi surface that is defined by a critical 
angle 9c measured at the (tt, tt) point from the Brillouin 
zone boundary towards the nodal directioni^ 



III. PENETRATION DEPTH 
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These two spectral functions enter the expressions for 
the penetration depth and the optical conductivity which 
are described in the following sections. In the limit of 
Age = but with a finite pseudogap, there are still two 
branches in the quasiparticle energy dispersion as shown 
in Fig. [TJb) for various symmetry directions in the first 
quadrant of the square Brillouin zone. The weighting fac- 
tors are shown as upright bars of varying height in 
the third dimension. Comparing with a case for no pseu- 
dogap where there is just the one energy dispersion with 
full weight [Fig.[IJa)], we see that a gap in energy due to 
the pseudogap formation has opened around the antin- 
odal direction in the Brillouin zone quadrant at (tt, 0) 
and (0, tt), and along the AFBZ from (tt, 0) to (0, tt). The 
presence of the two bands gives rise to interband tran- 
sitions in the frequency-dependent optical conductivity 
and a depletion of the intraband Drude component as we 
will later discuss. 

In the following work, we will also examine several 
other models in relation to the YRZ model. We will re- 



The London penetration depth A(r) is given by the 
zero frequency limit of the imaginary part of the optical 
conductivity (t{T, fi) at temperature T. Specifically, 



1 47r 



(7) 



where c is the velocity of light. The imaginary part 
of cr(T, n) has both a contribution from the paramag- 
netic and diamagnetic part of the current. A particularly 
transparent formal expression, which manifestly vanishes 
when the superconducting gap vanishes, results when 
both contributions are treated on the same footing^i^ 



1 



= lim 
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q^O \ Lu'' 
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where e is the electron charge and /(w) is the Fermi 
function 1/[1 -|- exp(/3a;)], with (3 = l/{kBT) and ks the 
Boltzmann constant. B(k,uj) is the anomalous spectral 
density previously introduced and is proportional to the 
superconducting gap. Consequently, this expression van- 
ishes in the normal state as it must. After some lengthy 
but straightforward algebra, one arrives at a more ex- 



4 



Weight 



Weight 



(0,0) 



(0,0) 



X = 0.20 




x = 0.14 




Ek/to 



pressed for now the factor of which weights this quan- 
tity in the YRZ model. We will be presenting all of our 
numerical work for \/\^{T) in units of e^to/H^d, where d 
is the c-axis distance per Cu-0 plane (not shown explic- 
itly in the supcrfluid density formulas given here which 
are written for 2 dimensions) . 

Our basic formula, Eq. can be reduced in sim- 
pler models such as the nodal liquid or the usual Fermi 
arc model. In both cases, the pseudogap is placed on 
the Fermi surface. For the arc model it is nonzero 
only on arcs about the antinodal direction. To place 
the pseudogap on the Fermi surface instead of having 
it on the antiferromagnetic Brillouin zone, we take 
to be equal to so that ~ ±Ek, and — 

(1 ± iklEk)/2, with Ek = Jil + A2 . It follows that 
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'^pg ' ^SC 



Es,W+ + W^=l, (h±)2 



(1 + Ek/Es)/2 = u', (v^y = (1 - Ek/Es)/2 = v', and 



^ uv ^ Asc/2Es- We then obtain 
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FIG. 1: (Color online) The electronic band structure for (a) 
X = 0.2 with no pseudogap and for (b) a:: = 0.14 with a pseu- 
dogap. In the YRZ model for the pseudogapped normal state, 
two bands, and , exist with fc-dependent weighting fac- 
tors and WjT , respectively, shown as a peak of height set 
by the weight. The EJ^ band is the lower energy band shown 
in heavy black and the E^ band is the higher energy band 
shown in lighter red. 



plicit formula for 1/X'^{T) in the YRZ model of the form 
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This formula is for the clean limit and ignores the in- 
elastic effects which are not large for the penetration 
depthf21 Only the coherent part of the Green's function 
contributes to the condensate. Note that we have sup- 



If we furth er take the pseudogap to be zero Es = 
+ then Eq. pI7|) becomes a well-known formula 
for the Fermi liquid (FL) penetration depth in BCS 
theory3^i2ii^i2£i^ With the pseudogap non-zero, we will 
refer to Eq. (|10p as the nodal liquid limit and if the pseu- 
dogap is cut off at a certain critical angle O,. away from 
the antinodal direction so that there is no pseudogap on 
an arc about the nodal direction, we will refer to this as 
the arc model. 

Examining various limits of Eq. (jlOp provides under- 
standing of the physics when a pseudogap is present. We 
recall that for the BCS case with no pseudogap, at zero 
temperature the tanh(/3i?/2) is equal to 1 and the T — 
limit of the penetration depth takes on a particularly 
simple form 
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— / de 



[A°,cos(2^?)]2 



{e2 Jr [A0,cos(26i)]2}3/2 



(11) 



where the continuum limit of the free electron bands has 
been taken to simplify the mathematics and A^(0) is the 
electronic density of states, assumed to be constant. The 
ratio of the electron density to the electron mass n/m — 
N{0)vp. Further, the last integral over energy in Eq. (fTTj) 
is independent of the factor of [A^^. cos(20)] and equal 
to 2. The angular integration normalized to 27r is then 
trivially equal to 1 . Thus we obtain the well-known result 
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that 1/A^(0) = Anne^ / mc^ . Moreover, in the hmit T 
0, the leading temperature dependence is given by the 
derivative term {—d/dE) in Eq. (fTO| . If we write for 
simphcity Ag = A^^ cos(20), we arrive at 



y{T) 'J, 2n 



de 



A^ 



(12) 



As T this last integral is strongly peaked about e ~ 
and As ~ (i.e., the nodal direction) and we obtain the 
standard result 
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For the case of the nodal liquid, the pseudogap is assumed 
to go like Apg(6') = A°g cos(26') over the entire Fermi sur- 
face as does the superconducting gap. In this case, ex- 
pression pUj) can be cast in the form of the standard BCS 
case with two changes. The square of the gap amplitude 
is to be replaced by the sum of the squares of supercon- 
ducting and pseudogap, i.e., A"^ — > ^ i^2c)^ + (^pg)^' 

and an overall factor of (A°JV[(^°c)^ + (^Pg)^] ^^"^ 
multiplies the entire expression. This leads immediately 
to the result 
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In this case the London penetration depth at T = is 
no longer simply a normal state property but depends 
explicitly on the value of the superconducting gap am- 
plitude. Also, the normalized slope of the linear-in- 
temperature term is greatly reduced when the pseudogap 
is large. While this limit is helpful because of its simplic- 
ity, we will see that it does not describe the YRZ results 
well. By contrast, the arc model which is closely related 
to the above equations does well in capturing the main 
physics contained in the more mathematically complex 
YRZ model. Before turning to this case, we note that at 
T = 0, when the pseudogap is small compared with the 
superconducting gap, the penetration depth is modified 
to 1/A2(0) = (47meV™c2)[l-(A0g/A3"j2] which tells us 
that the superfluid density is effectively reduced over its 
no pseudogap value. The pseudogap competes with the 
superconductivity for phase space. 

Introduction of an arc over which the pseudogap is 
zero, while finite from antinodal direction to 0c, mod- 
ifies Eq. pT|) but at the same time leaves the sec- 
ond integral in Eq. p2)) completely unaltered because. 



for sufficiently small temperature, only the angular re- 
gion very close to the nodal direction is of importance. 
But in these regions, the pseudogap is zero so that the 
first derivative term —d/dE in Eq. (fTU]) is unaffected 
by the pseudogap. Therefore, we obtain for this term 
{Anne'^ /mc'^)[—(2\n2)kBT / A'^^] which is completely un- 
changed from its Fermi liquid value. This is not so for the 
value of the zero temperature penetration depth. This 
quantity does know about the pseudogap. In the contin- 
uum approximation, it is given by 
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where for simplicity we have assumed the pseudogap to 
have the same angular dependence as the superconduct- 
ing gap, namely the d-wave cos(20) dependence. Here 
the bar on Apg is to mean that it is zero in the in- 
terval 9c to 7r/4 and all other symmetry related inter- 
vals. Next we note that the energy integral will give 
2(A0JV[(A0J2 + (A0g)2] from the regions where A^g is 
nonzero (antinodal) and 2 from the regions where A^g is 
zero (nodal). We finally get 
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For 6'c = (that is, no pseudogap), we get back the 
classical London result of Eq. pT|) after integration over 
energy and angle. For Oc = 7r/4, the nodal liquid result 
of Eq. (dH) in the T ^ limit is obtained. When the 
pseudogap A^g = 0, we recover the known result of Fermi 
liquid theory, that the zero temperature value of the pen- 
etration depth depends only on normal state parameters 
and not on the explicit value of the superconducting gap 
itself. When the pseudogap is small as compared to the 
superconducting gap, the correction to 1 in Eq. (|16p is 
small and of order [AOd n){A^^/ A^^)'^ . In the opposite 
limit of large pseudogap, the superconducting gap drops 
out and the correction is A9c/n, explicitly independent 
of both gaps. The pseudogap, however, implicitly deter- 
mines the critical angle dc related to the part of the Fermi 
surface over which the pseudogap is non-zero. 



IV. OPTICAL CONDUCTIVITY 

The real part of the optical conductivity, 'Rea{T, ft), is 
jiven by2» 
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In the clean limit of the YRZ model, we obtain, after long but straightforward algebra, the expression 
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The first term in Eq. (fT8|) is proportional to 5{Vl) and a 
Drude weight can be defined as 



^ Jo 



dE^ 



(19) 



The second piece is the interband contribution because 



it involves transitions between 



Et 



and Eg , and after 



integration over fl we have its weight Wjb equal to 
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(20) 



These two quantities are closely related to the pene- 
tration depth. As in that case, an interesting limit of 
these expressions occurs when the pseudogap is placed 
on the Fermi surface, i.e., taking — i^^, the same 
algebraic simplifications as previously described apply. 
First it is instructive to look at Eq. before inte- 

gration over photon energy Q. The second term van- 
ishes because the combination of u's and v's in the small 
square bracket vanishes and the third (last) term be- 
comes 2Wi^W^[—df{Eg)/dEg]S{fl) and combines with 
the first term to give 



dfjEg) 
dEs 



(21) 



where we have noted that + = 1. In the on-the- 
Fermi-surface limit of the pseudogap, the sum of the two 
contributions to the total weight can still be denoted by 



Wd/'2' as before and 



Wr 



dfjEg) 
dEg 



(22) 



In the continuum or free electron model for the band 
structure, with neither superconductivity nor pseudogap 
and at zero temperature, Wd = nne^ /m = rip/4, a well 
known result where f2p is the normal state plasma fre- 
quency. For superconductivity over the entire Fermi sur- 
face, but no pseudogap, we get zero as expected. In 
the clean limit at T = 0, the entire real part of the 
conductivity goes into the condensate. This would also 
hold for the nodal liquid. However, in the arc model 
with no superconductivity, we have a finite Wd, namely 
[n7re^/m](l— 46'c/7r) corresponding to the gapless arc that 
remains about the nodal direction. The limit 6c — Q gives 
no reduction from its Fermi liquid value and 6c = 7r/4 
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Doping X 

FIG. 2: (Color online) (a) The Fermi wavevector kp and ve- 
locity VF in the nodal direction for the YRZ dispersion. These 
are given in units of 1/a and toa/h, respectively, (b) The real 
part of the optical conductivity (arbitrary units) in the pseu- 
dogapped normal state as a function of fl/to, illustrating the 
Drude and interband components, (c) A plot comparing the 
plasma frequency Qp and the square of the inverse penetra- 
tion depth. Also, shown are the Drude and interband optical 
spectral weights, Wd and Wib as a function of doping x and 
the interband spectral weight in the superconducting state 
WiB,sc- For all of these quantities shown here we have left 
out the doping-dependent prefactor gt- These curves are all 
in units of e^to/h^d, where d is the c-axis distance per Cu-O 
plane. 



corresponds to a fully gapped Fermi surface (nodal liq- 
uid). 

In conventional BCS theory, the Ferrell-Glover- 
Tinkham^^i^ sum rule states that the optical spectral 
weight is not changed on entering the superconducting 
state. In terms of the notation introduced here, it reads 

2 poo pea 

-— — +/ dnReas{T,n)= dURea {T , il) , 

(23) 

where we denote by S/N the superconducting/normal 
state, respectively. We begin with a discussion of the 
validity (or lack thereof) of this rule in the YRZ model. 
We will work only at T = which is simplest and where 
the superconducting gap is fully developed and so has 
its maximum effect on charge dynamics. Thus we would 
expect that, if the presence of the pseudogap leads to a 
violation of Eq. it should be most noticeable here. 

In Fig. [2Ic), we plot several of the relevant quantities as 
a function of doping, x. The solid black dots give rip/8, 
i.e., the right hand entry of Eq. ([25)1 . The change in be- 
havior at a; = 0.2 corresponds to the QCP in our phase 
diagram and signals the emergence of a finite pseudogap. 
Since we are in the clean limit, and also in the normal 
state, rip/8 corresponds to Wd/2 (green dash-double- 
dotted curve) in the Fermi liquid region of the phase 
diagram x > 0.2. Below this doping, has two con- 



tributions Wd/2 from the Drude peak and Wjb (long- 
dashed blue curve) from the interband transitions. A plot 
which illustrates the two contributions to Kea{T = 0, il) 
is shown in Fig. [DJb). For this plot we used a formula 
for the conductivity which included elastic impurity scat- 
tering with quasiparticle scattering rate rj set equal to 
O.Olio- While the Drude and interband contributions are 
not completely separated, the two distinct contributions 
remain clearly identified. The Drude is centered about 
fi = and the interband piece is shifted to higher ener- 
gies with an onset just above f2 = 0.2to. The plot is for 
doping X = 0.14. This is a case where the pseudogap is 
larger than the superconducting gap with Apg/io = 0.18 
which corresponds roughly to the onset mentioned above. 
Returning to Fig. EJc), the sum of Wd/2 and Wib add 
up to ilp/8 in the normal state. As we have already 
mentioned, in the superconducting state at zero temper- 
ature in the clean limit, the entire Drude condenses into 
the superfluid and there is no contribution to the optical 
spectral weig ht J^ReasiT = 0,n)dn fr om this term, 
only the interband term remains. Its value as a func- 
tion of doping is denoted by Wib,sc and is given by the 
solid black curve. Comparing this with its normal state 
value (long-dashed blue curve), we see that at doping 
just below the QCP at a; = 0.2, this contribution is very 
small and so most of Wjb also goes into the condensate. 
However, as x decreases towards the more underdoped 
regime, much less of Wjb condenses and by a; ~ 0.1, this 
condensation has stopped. This behavior is expected and 
shows that the interband piece is less susceptible to con- 
densation into Cooper pairs than is the Drude and that 
this trend increases rapidly as the pseudogap energy be- 
comes large as compared with twice the superconducting 
gap energy. In our phase diagram, this boundary comes 
at about ^0.1 doping as we would expect. 

The dot-double-dashed purple curve gives 1/8A(0)^. 
When this is added to Wib,sCi we obtain, within our 
numerical accuracy, the i.e., the solid black circles. 

Thus we find that to this precision, the YRZ model also 
satisfies Eq. (j23p which is one of our important results 
about optical spectral weight distribution in addition to 
our observation that little of the interband contribution 
to the conductivity condenses into Cooper pairs when Ag^ 
becomes small as compared to the pseudogap Ap^. Fi- 
nally, in Fig.[5fa), we show our results for the Fermi mo- 
mentum in the nodal direction (solid black curve) as well 
as the corresponding Fermi velocity (dashed red curve). 
Our values agree well with those given in the work of 
YRZ and will be needed in a following section. 



V. RESULTS 

We will now present in detail our results and analysis 
for the penetration depth. Note that up until now, in 
our formulas and discussion, we have not included the 
Gutzwiller gt prefactor for weighting only the coherent 
part of the Green's function in Eq. ([1]). Due to the prod- 
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FIG. 3: (Color online) The temperature dependence of the 
penetration depth shown as 1/\^{T) versus T/to. Here, the 
prefactor, that arises from the definition of the coherent part 
of the Green's function, is included. Shown are curves for the 
Fermi liquid state (FL, solid black curve) and the modification 
which occurs when a pseudogap is present in the YRZ model 
(red dashed curve). The curve labeled YRZ mod. (green 
dash-dotted) is a modification of the YRZ model where the 
superconducting gap is only non-zero on the Fermi pocket. 
The three Fermi surfaces involved are illustrated schemati- 
cally in the inset (see text) including the critical angle 9^. 



uct of two spectral functions that enters the formulas 
for the penetration depth and the optical conductivity 
[Eq. ^ and Eq. HT]), respectively], 1/A2(T) and a{T, 9) 
will include a doping-dependent prefactor of . We now 
include it from here on forward as this factor reflects the 
highly correlated nature of the cuprates and is necessary 
for a proper comparison with experiment. 

In Fig. [31 we present our results for the temperature 
dependence of the penetration depth for a; = 0.12. Three 
cases are considered and are to be compared. The solid 
black curve labeled Fermi liquid is obtained when the 
pseudogap is set equal to zero and is for comparison with 
the short-dashed red curve which includes the pseudogap 
and is labeled YRZ. We note two striking differences be- 
tween these two sets of results. First, as a function of 
temperature, the FL curve shows concave down behav- 
ior. This differs slightly from the perhaps more classical 
curve for a c?-wave superconductor in a model, where the 
continuum approximation is used for the band structure 
and a cos(20) gap variation is assumed on the cylindrical 
Fermi surface, but is mainly due to our use of a gap ratio 
of 2A'^^/kBTc = 6 rather than the weak coupling BCS d- 
wave value of 4.28. The use of a larger value for this ratio 
is in keeping with experimental observation that finds it 
of order of 6 (and occasionally reported to be even larger). 
Our results also represent a generalization in which the 
superconducting gap varies over the entire Brillouin zone 
and the band structure is given in tight-binding with up 
to third nearest neighbor hopping. The same gap and 
bands are used for the red dashed curve which by con- 
trast shows slight concave upward variation as a function 



of temperature. Here a pseudogap is included and this 
has changed the usual large Fermi surface of Fermi liquid 
theory which is shown in one quadrant of the Brillouin 
zone in the top inset, to the Fermi surface shown in the 
middle inset. Here, the blue curve is the Luttinger pocket 
of the YRZ theory and the red line is its shadow exten- 
sion which represents a momentum contour of minimum 
approach to the "Fermi surface" in regions of momen- 
tum space where a gap exists so there are no true zero 
energy states. This reconstruction of the Fermi surface 
into a Luttinger pocket, which by its construction con- 
tains exactly x empty states (holes), has led to a large 
suppression of the zero temperature penetration depth as 
expected in our simplified Eq. The second striking 

feature of these two curves is that they have identical 
values of slope, as a function of temperature T, out of 
T = 0. The formation of the Luttinger pockets, as the 
pseudogap increases in the underdoped regime of x below 
the QCP dX X — 0.2 in our phase diagram, reduces the 
amount of Fermi surface that is available for pairing as in 
other competing order parameter scenarios^, but leaves 
the nodal region ungapped. The very low temperature 
excitations out of the ground state are confined to the 
cone very near zero energy, which exists in the nodal di- 
rection only, but these excitations do not sample directly 
the pseudogap and so the slope retains its Fermi liquid 
value. 

The final curve in Fig. [3] (green dash-dotted) is a case 
where we have cut off the superconducting gap outside 
the solid angle 1 — 46'c/7r which defines the region of the 
Luttinger pocket, with 9c shown in the lower inset. As 
we expect from the above arguments this does not change 
the low temperature slope but does affect the zero tem- 
perature value of the superfluid stiffness which is further 
reduced over its Fermi liquid value. Our motivation for 
applying this cutoff is the expectation that the supercon- 
ducting gap will form mainly on the part of the Fermi 
surface which remains ungapped. This idea is consistent 
with recent ARPES data'*^ and will be discussed further 
in a future paper—. 

Figure [3] further emphasizes the insensitivity of the low 
temperature slope of 1/X^{T) to the formation of the 
Luttinger pockets and any cut off in the superconduct- 
ing gap away from the nodal region. Shown in Fig. 2] is 
U)d/dT[\/\^{T)] vs T/Tc for the three cases iflustrated 
in Fig. [31 The line labels are the same. We see perfect 
agreement between the three curves below T/Tc ^ 0.3. 
With increasing temperature, the Fermi liquid case shows 
a downward trend while by contrast the two other non- 
Fermi liquid curves show the opposite. In the inset to 
the figure we show our results for the slope in the YRZ 
model (solid black dots) as a function of doping. Also 
shown for comparison is the formula 

dl TT VA 

where va = |VA,c(fc)|fc^ = {A°^jV2)\sm{kF^)\, with aU 
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FIG. 4: (Color online) Plot of slope of the penetration depth 
curves shown in Fig. [3] versus T/Tc- The inset shows a com- 
parison of the magnitude of the slope as T ^ with the 
formula ^■Kg'i[2\n2/'K\{vF /va)- This expression (solid line) is 
shown tested against the numerical data (dots) as a function 
of doping X. 
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FIG. 5: (Color online) The temperature dependence of the 
inverse square of the penetration depth for various dopings. 
The inset shows the values of the superconducting gap and 
pseudogap in units of to for varying doping parameter x. 



quantities defined in terms of the units stated in Fig. [2] 
and elsewhere. Here, the |sin(A:Fa;)| factor accounts for 
the fact that the superconducting gap amplitude sam- 
pled on the Luttinger pocket is somewhat smaller than 
the input value A^^ which enters the phase diagram and 
corresponds to the maximum gap in the Brillouin zone. 
Formula (|24|) is derived for a Fermi liquid with general 
band structure and we see here that it also applies to 
YRZ theory (Fig.SD. 

In Fig. El we show 1/\^{T) vs T/Iq for different val- 
ues of doping X. The phase diagram identifying the su- 
perconducting gap and pseudogap values used is shown 
for reference in the inset where we see the QCP at 
X — 0.2 corresponding to the onset of a finite pseudo- 
gap which progressively modifies the large Fermi sur- 
face of Fermi liquid theory to smaller Luttinger pock- 



ets. Several overall trends displayed in the 1/A^(T) 
curves shown in this figure are in qualitative agree- 
ment with two recent data sets, one for highly under- 
doped YBasCuaOe+y (YBCO)i^>ii and the other for the 
Bi-based cuprates (Bi:2212) Bi2.i5Sri.85CaCu208+5 and 
Bi2.iSri.9Cao.85Yo.i5Cu208+5 (Ref4^). In addition to 
the reduction in the T = value of 1/A^(r), there is 
also a trend for the linear temperature dependence out 
of T = to remain to higher values of reduced temper- 
ature T/Tf.. This contrasts with the concave downward 
behavior seen in overdoped and optimally doped cases^^ 
as Tc is approached. Such a trend towards linearity is 
clearly seen in the data of Huttema et al4^i^ who note 
a near linear dependence of the data almost all the way 
to Tc. However, they also find, in Fig. 3 of their paper, 
a turnaround to a law at low temperatures in agree- 
ment with the well-known effect of impurities in d-wave 
superconductors. While we have not included impurities 
in our work, we expect no modifications to this law to 
arise in the YRZ model and so agree with these results. 
However, because the YBCO data is in the deeply un- 
derdoped regime near the bottom of the superconducting 
dome, there could be additional effects that become im- 
portant such as fiuctuations,*^ not included here. Hence 
we turn to the data of Anukool et al4^ on Bi:2212 and 
in particular their Fig. 3(b) which, in agreement with 
our findings, shows near linear behavior over the entire 
temperature range considered. Similar behavior is also 
seen in their purer sample (shown in their Fig. 3(a)) al- 
though in that case the overdoped samples show a low 
temperature upturn which may indicate physics not in- 
cluded here. The Bi:2212 data set will be examined more 
closely in what follows as the doping range is more com- 
patible with the assumptions of the YRZ model. 

In Fig. [HI we show our results for the value of the zero 
temperature penetration depth as a function of doping 
X for three models. The solid black circles are based on 
a Fermi liquid model for the renormalized band which 
includes the Gutzwiller factors of the YRZ theory but 
no pseudogap. We see a smooth evolution with increas- 
ing superfluid density as doping increases. The solid red 
squares are for comparison and include the pseudogap 
which is finite below the QCP. A finite Apg leads to Fermi 
surface reconstruction as shown in the insets of Fig. [3] 
where the Luttinger pockets are seen in the middle and 
lower frame and are to be compared with the large Fermi 
liquid Fermi surface of the top panel. For x just below 
the QCP, the Fermi surface reconstruction can be even 
more complex than in Fig. [3] An example is shown in 
the inset of Fig. [Slfor x = 0.18, where we show Luttinger 
surface with the dashed curve being the AFBZ bound- 
ary. Note the pieces of occupied fc-space on the other 
side of the AFBZ boundary. In the main frame, the solid 
green diamonds include the pseudogap and, in addition, 
the superconducting gap is assumed to be non-zero only 
on the Luttinger pocket in an attempt to make the cal- 
culations more realistic. As we saw in Eq. (|16p with a 
simplified continuum model for the band structure, we 
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expect 1/A^(0) to drop as the Fermi surface arc is re- 
duced because the size of the Luttinger pocket shrinks. 
The sohd black curve corresponds to the product of the 
Fermi hquid value of 1/A^(0) with the factor 1 — 4:9c/ tt, 
the latter gives an approximate measure of the ratio of 
the angle of the remaining Fermi arc compared with that 
for the Fermi liquid. This curve follows the same gen- 
eral trend as do the results of detailed calculations. This 
shows that Fermi surface reconstruction is responsible for 
the drop in solid green diamonds below the Fermi liquid 
values (solid black circles). We note again in this regard 
that as the QCP is approached from below {x < 0.2), 
there is a region where the Fermi surface reconstruction 
is complex (see insert) and is not easily characterized by 
a single arc length. Here, we use the values previously 
derived within the context of an application of the YRZ 
model to the specific heat^^. This explains why the solid 
black line starts to deviate from the solid green diamonds 
in the region near the QCP. We emphasize that the solid 
green diamonds include a cut off on the superconducting 
gap. 

As can be seen in the approximate qualitative formula 
(fTH]) . we expect the superconducting gap to affect the 
zero temperature condensate, a result which is quite dif- 
ferent from the Fermi liquid case where 1/A^(0) depends 
only on normal state parameters and not on the gap. Re- 
cently Kopnin and Sonin^ found a similar dependence 
on superconducting gap in the case of graphene although 
no superconductivity has yet been reported in this sys- 
tem even though carbon nanotubes have been found to 
super conduct 4^ Returning to Eq. the correction 

factor of c = (-46l,/7r)(A0g) V[(A°c)^ + i^pg?] reduces 
to —4:9c/tt only when A^^, is assumed to be zero when Apg 
is finite, the case described so far. If the superconducting 
gap is allowed to exist over the entire Fermi surface, the 
correction factor c will be less than —A9c/tt and we get the 
blue dash-dotted curve in Fig. [S] as the product of (1 — c) 



times the Fermi liquid value of 1/A(0)^ 



Eq. nil)]. 



This curve agrees well with the solid red squares which 
are obtained from the full calculations. It is clear that 
Fermi surface reconstruction leads to a strong reduction 
of the zero temperature superfluid density as a result of 
the loss of ungapped states on the Fermi surface. 

We now turn to experiment to assess the validity of 
the theoretical model due to YRZ. Recent results as a 
function of doping on the penetration depth of Bi-based 
cuprates (Bi:2212) have been published by Anukool et 
al.'*^. They give results for a range of doping from un- 
derdoped to overdopcd and for Tc's reduced by up to 
50%. Because their optimal doping is at p = 0.16 and 
our model uses x = 0.2, we shifted the data set doping 
values by 0.04 in order to match their dome with ours. 
Normalizing to the optimal doping, we then find that 
their curve matches our Tc dome if we convert our en- 
ergy gap dome to one for Tc by using a value of the gap 
ratio 2A'^^/kBTc = 6 and then taking to = 165 meV (our 
energy units have all been in units of to). With this we 
find a good match to the experimental Tc versus doping 
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FIG. 6: (Color online) The T — inverse square penetration 
depth as a function of doping, comparing the cases of the 
Fermi liquid (FL), YRZ and the modified YRZ ansatz. The 
shade region emphasizes the reduction in this quantity due 
to Fermi surface reconstruction. The curves are for simple 
analytic formulas as discussed in the text. The inset shows 
the Luttinger area (shaded blue) for the case of a; = 0.18, il- 
lustrating an intermediate regime of the Fermi surface recon- 
struction as a function of doping. This inset also illustrates 
the critical angle Oc for defining the Fermi pocket, which has 
become ill-defined in this case. 
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FIG. 7: (Color online) Tc vs doping and vs the superfiuid 
density normalized to the optimal value in the YRZ model. 
The solid blue curve is the result of the YRZ model and the 
red dashed curve is the modified YRZ model as discussed 
in the text. Comparison has been made with the data of 
Anukool et al.^^ with the Bi:2212 (solid dots) and Bi(Y):2212 
(open dots). 



curve shown in Fig. [71 This now fixes our parameters for 
the rest of our work. In Fig. [71 (right hand frame) we 
show the classic Uemura plot of Tc versus 1/A^(0), where 
we have normalized our results and the data to the value 
at optimal doping. The agreement with the data is ex- 
cellent and on the undcrdoped side, the data appears to 
favor more the modified YRZ calculation. Furthermore, 
in Fig. [51 we show the data compared with the inverse 
penetration depth squared versus doping and once again 
the data favors the modified YRZ curve. What is clear 
here is that the data agrees well with the Fermi liquid 
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FIG. 8: (Color online) Comparison of the superfluid density 
at T = with the Bi:2212 (solid dots) and Bi(Y):2212 (open 
dots) data of Anukool et al.— . This is shown as a function 
of doping scaled to the optimal doping value. The solid blue 
curve is the result of the YRZ model and the red dashed curve 
is the modified YRZ model as discussed in the text. The inset 
shows the slopes of the data as T — > in comparison with the 
analytic formula discussed in Fig. U 



curve above optimal doping and on the underdoped side 
it follows the YRZ curve which includes the reconstructed 
Fermi surface. Indeed, the drop in the data below opti- 
mal doping can be taken as possible indirect experimen- 
tal evidence for Fermi surface reconstruction as the Fermi 
liquid curve would have been higher (as discussed in re- 
lation to Fig. [5]). All is not ideal, however. We have used 
a normalized quantity here. In absolute value, we find 
that the superfluid density is off by about a factor of 3- 
4. Nonetheless, this can be rectified as the penetration 
depth depends only on to and so one could increase the 
value of to to account for this. This would then require an 
adjustment of the gap ratio to keep the good agreement 
with the Tc dome of Fig. [3 but we find that it would be 
unrealistically large. Another possibility is to note that 
our estimate of the Fermi velocity in the YRZ model at 
optimum is low by a factor of two. The superfluid den- 
sity shown here goes as the square of the velocity and 
this could possibly correct the situation. A related issue 
is that the band structure used in the YRZ paper is for 
a very different compound than that of Bi:2212 and so 
this could change some of these quantitative numbers. 
We wish to stress that the merit of this theory should be 
seen in its ability to give qualitative insight into the pseu- 
dogap phase and the good agreement that we find here 
is very encouraging given the lack of detailed parameter 
fitting to this particular material. Also shown in Fig. [8] 
is an inset plotting the experimental slope of the inverse 
square penetration depth curves for T —f 0. Again this 
has been normalized to the optimal doping value and we 
find good agreement with the analytic formula discussed 
for the slope in Fig. U) An analysis based on ARPES has 
given somewhat different results4^ It should be noted 
that we extracted the slope values ourselves and so this 




FIG. 9: (Color online) A comparison of the effect that various 
models have on the square of the inverse penetration depth as 
a function of temperature. The schematic diagrams illustrate 
the momentum dependence of the two energy gaps for the 
Fermi arc model and its modification, relative to the nodal 
liquid case. 



is not a rigorous representation of the data and indeed 
the temperature-dependence of some of the experimental 
curves was unusual in a few cases, but we took the lowest 
temperature value of the slope in any event. The optimal 
doping curve was in this category as it showed an upturn 
at low temperature in the case of the pure sample and so 
its slope value and a few of the others could be revised in 
the hands of experimental experts. However, the inset of 
Fig. [8] (as well as the other comparisons of Figs. [7] and [8|) 
serves to illustrate that the Gutzwiller prefactor for the 
coherent part of the Green's function, leading to a pref- 
actor of , appears to be essential in giving the correct 
trend of the data with doping. This points to very strong 
correlations in these systems. 

In Fig[9l we give results which show the relationship of 
the YRZ model with two other prominent models. Both 
involve the assumption that the pseudogap acts on the 
Fermi surface. This corresponds to taking the limit of 
~^ Cfc i'^ '^^^ Elq. (HI) for the penetration depth, i.e., 
to replacing the AFBZ energy by the carrier dispersion 
curve. As we saw, this greatly simplifies the equation and 
reduces it to the familiar form of Eq. ^TU\i for a BCS su- 
perconductor with one modification. While in Eq. 
the superfluid density remains directly proportional to 
a factor of A^^ and so manifestly vanishes in the nor- 
mal state, the energy Es = \J^k + ^sc + ^pg involves 
the sum of the squares of Age and Apg rather than the 
square of Age alone. It is only through this factor that 
the pseudogap enters in this simplification. The Fermi 
arc model^S. corresponds to taking Apg non-zero only on 
an arc centered about the antinodal direction, leaving 
an ungapped region centered around the nodal direction 
(on the Luttingcr pocket) as illustrated by the center in- 
sert in Fig. [51 One can construct realistic modelsi^ for 
the momentum variation of the pseudogap on the Fermi 
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surface. For simplicity, here we retain its cosine form 

A" {x) 

Apg = — (cosfca;a — COS kyu) but cut it off on an arc 
defined by the angle 0c previously introduced. Taking 
dc,aic — 24°, chosen to get a best fit to the penetration 
depth curve of the YRZ model (solid black curve), one 
obtains the red open circles which overlap the black solid 
curve almost perfectly. Thus with the right choice of 
^c.arc, one can get an almost perfect match between the 
two models. If we further assume that the superconduct- 
ing gap is present only on the ungapped (by the pseu- 
dogap that is, see right hand inset) part of the Fermi 
surface, we obtain the dashed curve denoted Arc mod. 
This has pushed the penetration depth down by about 
the same amount as we saw in Fig. |3] for the full YRZ cal- 
culation when we included a superconducting gap mod- 
ification in this model. The choice of 9c. arc in Fig. [9] is 
important as can be seen from the dash-dotted blue curve 
which is also based on an arc model but with a different 
cut off model 9c,yrz = 26°, a value obtained from the 
size of the arc subtended by the actual Luttinger pocket 
of the YRZ theory. While this choice decreases further 
the magnitude of the superfluid density because there is 
even less ungapped arc, it does not change the qualitative 
behavior obtained. 

Finally, the solid green curve with triangles is our re- 
sult for the nodal liquid which corresponds to lifting the 
cut off on the pseudogap. In this case an effective gap 

of ^ A^^ -I- A^g replaces the usual superconducting gap 

in the standard Eq. pO)) for 1/A^(r) with one critical 
difference. In the numerator of Eq. (fTUI) . it is still A^^ 
which remains and not A^„-|-A^„. Since we have assumed 
the same momentum dependence for the superconduct- 
ing gap and pseudogap, we can replace the explicit Ag^ 
factor in Eq. (fTO|) by A^^ -I- Apg and take out the fac- 
tor of (A°J2/[(A0J2 _^ (A°g)2] to compensate for this. 
Thus the equation for the penetration depth becomes 
the standard one for an effective superconducting gap of 

y^A^~+~A^ but with the difference that a factor of the 

square of superconducting-to-effective-gap amplitude is 
to multiply the entire expression as we see explicitly in 
our simplified expression (|14p . While these simplifica- 
tions allow us to obtain simple analytic expressions, we 
see that this limit fails to give quantitative results when 
compared with YRZ. Compare the green curve with tri- 
angles with the solid black curve in Fig [9] In fact one 
should not expect the two models to agree. The nodal 
liquid is the limit when the size of the Luttinger pocket 
tends to zero. But this is never reached inside the super- 
conducting dome. In our calculations at a doping where 
Tc has been depressed to zero, there remains a sizable 
part of the Fermi surface which is ungapped. It is not 
surprising then that the nodal liquid should show quali- 
tative behavior not part of the arc model, such as, a su- 
perfluid density at T = which decreases by a factor of 
the square of the ratio of superconducting to pseudogap. 
The slope of the linear-in-T variation at low tempera- 



ture also becomes inversely proportional to y A^g -|- A^^ 
rather than to Age and so becomes very flat. 

VI. CONCLUSIONS 

YRZ have provided a simple model for the coherent 
part of the charge carrier Green's function which applies 
below a quantum critical point characterizing the be- 
ginning of pseudogap formation. As the Mott-Hubbard 
transition to an insulating state is approached with de- 
creasing doping, the magnitude of the pseudogap in- 
creases, the bands narrow and the weight of the coher- 
ent piece of the Green's function decreases according to 
well-defined Gutzwiller factors which account for correla- 
tion effects. With increasing pseudogap magnitude, the 
Fermi surface reconstructs. It goes from the large Fermi 
surface of Fermi liquid theory to ever smaller Luttinger 
pockets. This leads directly to a reduction of the value of 
the zero temperature inverse squared penetration depth 
because there are fewer ungapped states which are avail- 
able to form the condensate. We show that this reduc- 
tion is roughly proportional to the ratio of the remaining 
arc length, defined by the Luttinger pocket, to the full 
length of the corresponding Fermi liquid surface. On the 
other hand, and in sharp contrast, the coefficient of the 
linear-in-temperature term of 1 /A^(T) at small T remains 
largely unaffected because this quantity depends only on 
the available very low energy excitations and these are 
confined to the vicinity of the nodal Dirac points. This 
region is not importantly affected by pseudogap forma- 
tion and implied Fermi surface reconstruction. 

A comparison of our results with recent experimental 
data on Bi:2212 gives good qualitative agreement and 
demonstrates the importance of the strong dependence 
on doping of the coherence weight in the YRZ model 
which derives from strong correlation effects. While this 
is a simple multiplicative factor, it accounts for a signifi- 
cant part of the reduction in superfluid density as the end 
of the superconducting dome is approached in the deeply 
underdoped region of the phase diagram. An additional 
reduction is due to Luttinger pocket formation which 
starts at the doping associated with the quantum crit- 
ical point and this effect provides a more abrupt change 
in magnitude which should be measurable as a signature 
of Fermi surface reconstruction. 

The well-known Ferrell-Glover-Tinkham sum rule of 
conventional superconductivity was found to apply here 
as well. The optical weight lost on entering the super- 
conducting state below Tc reappears in its entirety in the 
superfluid fraction. 

Comparison of our results with those based on an 
arc model, with pseudogap placed on the Fermi surface 
rather than on the AFBZ of the YRZ model, can be made 
to agree very well if the arc length on which the pseu- 
dogap is assumed to be non-zero is appropriately cho- 
sen to obtain a best fit with the YRZ case. The nodal 
liquid concept corresponds to the limit when the Fermi 
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surface is fully gapped by the pseudogap except for the 
nodal points. This limit is never reached in YRZ theory 
because the size of the Luttinger pocket remains quite 
significant at the doping which corresponds to the end 
of the superconducting dome. Nevertheless, the nodal 
liquid limit remains valuable because it yields analytic 
results which can provide useful insight into the deeply 
underdoped case. 



Acknowledgments 



This work has been supported by NSERC of Canada 
and by the Canadian Institute for Advanced Research 
(CIFAR). 



* Electronic address: |nicol@physics.uoguelph.ca| 

^ T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). 
^ See J. P. Carbotte, E. Schachinger, and D. N. Basov, Na- 
ture (London) 401, 354 (1999) and references therein. 
^ E. Schachinger and J. P. Carbotte, Phys. Rev. B 62, 9054 

(2000) . 

* S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. Nayak, 
Phys. Rev. B 63, 094503 (2001). 

^ V. J. Emery and S. A. Kivelson, Nature (London) 374, 
434 (1995). 

^ M. Randeria, N. Trivedi, A. Moreo, and R. T. Scalettar, 
Phys. Rev. Lett. 69, 2001 (1992). 

G. Alvarez and E. Dagotto, Phys. Rev. Lett. 101, 177001 
(208). 

^ Q. Chen, K. Levin, and I. Kosztin, Phys. Rev. B 63, 184519 

(2001) . 

^ S. Hiifner, M. A. Hossain, A. Damascelli, and G. A. 

Sawatzky, Rep. Prog. Phys. 71, 062501 (2008). 
^° H. V. Lohneysen, A. Rosch, M. Vojta, and P. Wolfe. Rev. 
Mod. Phys. 79, 1015 (2007). 

L. Taillefer, arXiv:0901.2313 and references therein. 

For example, A. Kanigel, U. Chatterjee, M. Randeria, M. 

R. Norman, S. Souma, M. Shi, Z. Z. Li, H. Raffy, and J. 

C. Campuzano, Phys. Rev. Lett. 99, 157001 (2007); A. 
Kanigel, M. R. Norman, M. Randeria, U. Chatterjee, S. 
Souma, A. Kaminski, H. M. Fetwell, S. Rosenkranz, M. 
Shi,T. Sato, T. Takahashi, Z. Li, H. Raffy, K. Kadowaki, 

D. Hinks, L. Ozyuzer, and J. C. Campuzano, Nat. Phys. 
2, 447 (2006); T. Kondo, R. Khasanov, T. Takluchi, J. 
Schmalian, and A. Kaminski, Nature (London) 457, 296 
(2009). 

" J. Meng, G. Liu, W. Zhang, L. Zhao, H. Liu, X. Jia, D. 
Mu, S. Liu, X. Dong, W. Lu, G. Wang, Y. Zhou, Y. Zhu, 
X. Wang, Z. Xu, C. Chen, and X. J. Zhou, arXiv:0906. 26 82. 
J. Hwang, J. P. Carbotte, and T. Timusk, EPL 82, 27002 
(2008). 

J. G. Storey, J. L. Tallon, G. V. M. Williams, and J. W. 
Loram, Phys. Rev. B 76, 060502(R) (2007). 
Y. Kohsaka, C. Taylor, P. Wahl, A. Schmidt, Jhinhwan 
Lee, K. Fujita, J. W. AUdredge, Jinho Lee, K. McElroy, 

H. Esaki, S. Uchida, D.-H. Lee, and J. C. Davis, Nature 
(London) 454, 1072 (2008). 

K.-Y. Yang, T. M. Rice and F. C. Zhang Phys. Rev. B 73, 
174501 (2006). 

P. W. Anderson, Science 235, 1196 (1987). 

B. Valenzuela and E. Bascones, Phys. Rev. Lett. 98, 

227002 (2007). 

^° E. lUes, E. J. Nicol, and J. P. Carbotte, Phys. Rev. B 79, 
100505(R) (2009). 

K.-Y. Yang, H. B. Yang, P. D. Johnson, T. M. Rice and 
F. C. Zhang, EPL 86, 37002 (2009). 



J. P. F. LeBlanc, E. J. Nicol, and J. P. Carbotte, Phys. 
Rev. B 80, 060505(R) (2009). 

W. N. Hardy, D. A. Bonn, D. C. Morgan, R. Liang, and 

K. Zhang, Phys. Rev. Lett. 70, 3999 (1993). 
'^^ Y. Uemura et al., Phys. Rev. Lett. 62, 2317 (1989). 
25 P. A. Lee and X.-G. Wen, Phys. Rev. Lett. 78, 4111 (1997). 

D. E. Sheehy, T. P. Davis, and M. Franz, Phys. Rev. B 70, 

054510 (2004). 

2^ D. Branch and J. P. Carbotte, Phys. Rev. B 52, 603 (1995). 
2^ P. G. Tomlinson and J. P. Carbotte, Phys. Rev. B 13, 4738 
(1976). 

2^ H. K. Leung, J. P. Carbotte, D. W. Taylor, and C. R. 
Leavens, Can. J. Phys. 54, 1585 (1976). 
H. K. Leung, J. P. Carbotte, and C. R. Leavens, J. Low 
Temp. Phys. 24, 25 (1976). 

D. Branch, Ph.D. thesis, McMaster University, 1997 (un- 
published) 

32 R. J. Radtke, V. N. Kostur, and K. Levin, Phys. Rev. B 
53, R522 (1996). 

33 E. Schachinger, J. P. Carbotte, and F. Marsiglio, Phys. 
Rev. B 56, 2738 (1997). 

^'^ C. O'Donovan and J. P. Carbotte, Phys. Rev. B 52, 16208 
(1995). 

C. O'Donovan and J. P. Carbotte, Phys. Rev. B 52, 4568 

(1995) . 

C. O'Donovan and J. P. Carbotte, Physica C 252, 9433 

(1996) . 

3^ S. Misawa, Phys. Rev. B 51, 11791 (1995). 

3** R. A. Ferrell and R. E. Glover, Phys. Rev. 109, 1398 

(1958) . 

3^ M. Tinkham and R. A. Ferrell, Phys. Rev. Lett. 2, 331 

(1959) . 

*° S. G. Sharapov and J. P. Carbotte, Phys. Rev. B 73, 
094519 (2006). 

T. Kondo, R. Khasanov, T. Takeuchi, J. Schmalian, and 
A. Kaminski, Nature (London) 457, 296 (2009). 
"2 J. P. F. LeBlanc, J. P. Carbotte, and E. J. Nicol, (unpub- 
lished) . 

■'3 D. M. Broun, W. A. Huttema, P. J. Turner, S. Ozcan, B. 
Morgan, R. Liang, W. N. Hardy, and D. A. Bonn, Phys. 
Rev. Lett. 99, 237003 (2007). 

W. A. Huttema, J. S. Bobowski, P. J. Turner, R. Liang, 
W. N. Hardy, D. A. Bonn, and D. M. Broun, Phys. Rev. 
B 80, XXX (2009). 

W. Anukool, S. Barakat, C. Panagopoulos, and J. R. 
Cooper, Phys. Rev. B 80, 024516 (2009). 
'^^ I. Hetel, T. R. Lemberger, and M. Randeria, Nature Phys. 
3, 700 (2007). 

■'^ N. B. Kopnin and E. B. Sonin, Phys. Rev. Lett. 100, 
246808 (2008). 

Z. K. Tang, Lingyun Zhang, N. Wang, X. X. Zhang, G. H. 



14 



Wen, G. D. Li, J. N. Wang, C. T. Chan, and Ping Sheng, 
Science 292, 2462 (2001). 

J. Mesot, M. R. Norman, H. Ding, M. Randeria, J. C. Cam- 
puzano, A. Paramekanti, H. M. Fretwell, A. Kaminski, T. 



Takeuchi, T. Yokoya, T. Sato, T. Takahashi, T. Mochiku, 
and K. Kadowaki, Phys. Rev. Lett. 83, 840 (1999). 



